Formation of ultracold SrYb molecules in an optical lattice by pho- 
toassociation spectroscopy: theoretical prospects 
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State-of-the-art ab initio techniques have been appUed to compute the potential energy curves for the SrYb molecule in the 
Bom-Oppenheimer approximation for the ground state and first fifteen excited singlet and triplet states. All the excited state 
Q_i potential energy curves were computed using the equation of motion approach within the coupled-cluster singles and doubles 
<^ framework and large basis-sets, while the ground state potential was computed using the coupled cluster method with single, 
double, and noniterative triple excitations. The leading long-range coefficients describing the dispersion interactions at large 
interatomic distances are also reported. The electric transition dipole moments have been obtained as the first residue of the 
polarization propagator computed with the linear response coupled-cluster method restricted to single and double excitations. 
I— I Spin-orbit coupling matrix elements have been evaluated using the multireference configuration interaction method restricted to 



Oh' 



single and double excitations with a large active space. The electronic structure data was employed to investigate the possibility 



I of forming deeply bound ultracold SrYb molecules in an optical lattice in a photoassociation experiment using continuous-wave 

^ lasers. Photoassociation near the intercombination line transition of atomic strontium into the vibrational levels of the strongly 

Q spin-orbit mixed b^T.^, a^Yl, A'n, and C'n states with subsequent efficient stabilization into the v" — 1 vibrational level of 

+^ 'the electronic ground state is proposed. Ground state SrYb molecules can be accumulated by making use of collisional decay 

• from v" = 1 to y" = 0. Alternatively, photoassociation and stabilization to v" = can proceed via stimulated Raman adiabatic 

Q passage provided that the trapping frequency of the optical lattice is large enough and phase coherence between the pulses can 

• 'tJ 'be maintained over at least tens of microseconds. 

^ 1 Introduction 

Molecules cooled to temperatures below T = 10^^ K allow 
for tackling questions touching upon the very fundamentals 
t> of quantum mechanics. They are also promising candidates 
in novel applications, ranging from ultracold chemistry and 
precision measurements to quantum computing. Cold and ul- 



o 



tracold molecules are thus opening up new and exciting areas 
of research in chemistry and physics. Due to their permanent 
electric dipole moment, polar molecules are particularly in- 
teresting objects of study: Dipole-dipole interactions are long 
7— I range and can precisely be controlled with external electric 
fields turning the experimental parameters field strength and 
• ^ orientation into the knobs that control the quantum dynam- 
ics of these molecules. It is not surprising then that a major 
objective for present day experiments on cold molecules is to 
achieve quantum degeneracy for polar molecules. Two ap- 
proaches to this problem are used - indirect methods, in which 
molecules are formed from pre-cooled atomic gases, '"^ and 
direct methods, in which molecules are cooled from molecu- 
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lar beam temperatures, typically tens of Kelvin.—— 

Direct cooling techniques, based on buffer gas cooling— 
or Stark deceleration^, produce cold molecules with a tem- 
perature well below 1 K. However, a second-stage cooling 
process is required to reach temperatures below 10^^ K. The 
most promising second-stage technique is sympathetic cool- 
ing where cold molecules are introduced into an ultracold 
atomic gas and equilibrate with it. Sympathetic cooling has 
successfully been used to achieve Fermi degeneracy in ^Li— 
and Bose-Einstein condensation in "^'K^' and to obtain ultra- 
cold ions. '^"'^ For molecular systems, however, sympathetic 
cooling has not yet been attempted, and there are certainly 
many challenges to overcome. In fact, calculations of the scat- 
tering cross sections for the collisions of molecules that can be 
decelerated with ultracold coolant atoms suggest that sympa- 
thetic cooling will not be efficient in many cases.— 

Alternatively, direct methods first cool atoms to ultralow 
temperatures and then employ photoassociation^ or mag- 
netoassociation^ to create molecules, reaching translational 
temperatures of the order of a few /iK or nK. In particu- 
larly fortuitious cases, photoassociation may directly produce 
molecules in their vibrational ground state.ii. Typically, how- 
ever, the molecules are created in extremely weakly bound 
levels, and follow-up stabilization to the ground state is nec- 
essary. For molecules built of alkali metal atoms, this has 
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been achieved using stimulated emission pumping ■ or alter- 
natively, employing coherent control techniques such as Stim- 
ulated Raman Adiabatic Passage (STIRAP^i^"— and vibra- 
tional cooling of molecules with amplitude-shaped broadband 
laser light. ^"^ 

Closed-shell atoms such as alkali earth metals are more 
challenging to cool and trap than open-shell atoms such as the 
alkalis. Closed-shell atoms do not have a magnetic moment 
in their ground state that enables magnetic trapping. More- 
over, for alkaline earth metals the short lifetime of the first ex- 
cited 'f state implies rather high Doppler temperatures, mak- 
ing dual-stage cooling a necessity where the second stage op- 
erates near an intercombination line. Despite these obstacles, 
cooling of calcium, strontium, and ytterbium atoms to micro- 
Kelvin temperatures has been achieved, and Bose-Einstein 
condensates of ^OCa,^^ 84^^,26^ 86^^,28^ 88^^,29 i70Yb,iO 

'^^Ybll have been obtained. 

In contrast to alkali metal dimers— or molecules consisting 
of an alkali metal atom plus a closed-shell atom,— the mag- 
netoassociation of two closed-shell atoms is not feasible ex- 
perimentally even if the nuclear spin is non-zero. The zero- 
field splittings and couplings between the atomic threshold 
and molecular states provided by the largest non-zero terms 
in the fine structure and hyperfine structure Hamiltonian for 
the electronic ground state, i.e., the scalar and tensor inter- 
actions between the nuclear magnetic dipole moments are 
simply too small. ^ On the other hand, the closed-shell struc- 
ture of the alkali earth metal and ytterbium atoms leads to 
very simple molecular potentials with low radiative losses and 
weak coupling to the environment. This opens new areas of 
possible applications, such as manipulation of the scattering 
properties with low-loss optical Feshbach resonances, high- 
resolution photoassociation spectroscopy at the intercombina- 
tion line, ^^'^^ precision measurements to test for a time vari- 
ation of the proton-to-electron mass ratio, quantum com- 
putation with trapped polar molecules, and ultracold chem- 
istry—. 

To the best of our knowledge, production of heteronuclear 
diatomic molecules built of closed-shell atoms has not yet 
been achieved experimentally. Also such processes have not 
yet been considered theoretically. Here we fill this gap and 
report a theoretical study of the photoassociative formation 
of heteronuclear diatomic molecules from two closed-shell 
atoms on the example of the SrYb molecule. Although the 
SrYb molecule may seem very exotic, especially for chemists, 
strontium and ytterbium atoms are promising candidates for 
producing molecules since they have both successfully been 
cooled and trapped. Moreover, both Sr and Yb have many 
stable isotopes. Such a diversity of stable isotopes is key 
to controlling the collisional properties of bosonic molecules 
with no magnetic moments and hyperfine structure. For ex- 
ample, one can effectively tune the interatomic interactions 



by choosing the most suitable isotope to achieve scattering 
lengths appropriate for evaporative cooling. Last but not least, 
we consider photoassociative formation of SrYb molecules 
since there are ongoing experiments'^" aiming at producing 
cold SrYb molecules in an optical lattice. 

The plan of our paper is as follows. Section |2]describes the 
theoretical methods used in the ab initio calculations and dis- 
cusses the electronic structure of SrYb in terms of the ground 
and excited-state potentials, transition moments, spin-orbit 
and nonadiabatic angular couplings. Section [3] analyzes the 
vibrational structure of the SrYb molecule as a prerequisite to 
determine an efficient route for photoassociation followed by 
stabilization into the vibronic ground state. It also discusses 
prospects of producing cold SrYb molecules by photoassoci- 
ation near the intercombination line of strontium, and subse- 
quent spontaneous or stimulated emission. Section |4] summa- 
rizes our findings. 

2 Electronic structure of SrYb 

In the present study we adopt the computational scheme suc- 
cessfully applied before to the ground and excited states of 
the calcium dimer— and of the (BaRb)+ molecular ion.— 
The potential energy curves for the ground and excited states 
of the SrYb molecule have been obtained by a supermolecule 
method: 

where denotes the energy of the dimer computed using 
the supermolecule method SM, and E^, X=A or B, is the 
energy of the atom X. For the ground state potential we used 
the coupled cluster method restricted to single, double, and 
noniterative triple excitations, CCSD(T). Calculations on all 
the excited states employed the linear response theory (equa- 
tion of motion) within the coupled-cluster singles and doubles 
(LRCCSD) framework The CCSD(T) and LRCCSD calcu- 
lations were performed with the DALTON program'*^. Note 
that the methods used in our calculations are strictly size- 
consistent, so they ensure a proper dissociation of the elec- 
tronic states, and a proper long-range asymptotics of the corre- 
sponding potential energy curves. This is especially important 
when dealing with collisions at ultralow temperatures, since 
the accuracy of the potential in the long range is crucial. 

For each electronic state we have also computed the long- 
range coefficients describing the dispersion interactions from 
the standard expressions (see, e.g. Refs. kslRgh that can 
be derived from the multipole expansion of the interatomic 
interaction operator The long-range dispersion coefficients 
were computed with the recently introduced explicitly con- 
nected representation of the expectation value and polariza- 
tion propagator within the coupled cluster method ^"-^^ and 
the best approximation XCCSD4 proposed by Korona and 
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Fig. 1 Potential energy curve (upper panel) and permanent dipole 
moment (lower panel) of the electronic ground state of the 

SrYb molecule. 

collaborators 52,. For the singlet and triplet states dissociat- 
ing into Sr(^P)H-Yb(i5), and Sr(^P)H-Yb('5), respectively, the 
dispersion coefficients were obtained from the sum-over-state 
expression with the transition moments and excitation ener- 
gies computed with the multireference configuration interac- 
tion method limited to single and double excitations (MRCI). 
The transitions from the ground X'E+ state to the and 
states are electric dipole allowed. The transition dipole 
moments for the electric, t/,, transitions were computed from 
the following expression—: 

rf,--(x'E+|d|(«)'|A|), (2) 

where d denotes the dipole moment operator Note that in Eq. 
^i^xoiy coiTesponds to transitions to ' n states, while / = z 
coiTesponds to transitions to 'E+ states. In the present calcu- 
lations the electiic transition dipole moments were computed 
as the first residue of the LRCCSD linear response function 
with two electiic, r, operators'*^. In these calculations we have 
used the DALTON program'*^. We have evaluated the depen- 
dence of the transition dipole moments with the internuclear 
distance for the same set of distances as the excited state po- 
tential energy curves. 

The electronic states of the low lying excited states of the 
SrYb molecule are coupled by nonadiabatic couplings. There- 
fore, in this work we have computed the most important angu- 
lar coupling matrix elements: 

A(/;) = ((«)2^+'|A||L+|(«')'^+'|A'|). (3) 

In the above expression denotes the raising electronic an- 
gular momentum operator Note that the electronic angular 
momentum operator couples states with A differing by one. 



9 [- -| 9 [- ^ \.'a'n - 

6 8 10 12 14 16 6 8 10 12 14 16 

R (bohr) R (bohr) 

Fig. 2 Potential energy curves of singlet (left) and triplet (right) 
excited states of SrYb dimer. 



Nonadiabatic couplings were obtained with the MRCI method 
and the MOLPRO code^^. We have evaluated the dependence 
of the nonadiabatic coupling matrix elements with the internu- 
clear distance for the same set of distances as the excited state 
potential energy curves. 

Strontium and ytterbium are heavy atoms, so the electronic 
states of the SrYb molecule are strongly mixed by the spin- 
orbit (SO) interactions. Therefore, in any analysis of the for- 
mation of the SrYb molecules the SO coupling and its depen- 
dence on the internuclear distance R must be taken into ac- 
count. We have evaluated the spin-orbit coupling matrix el- 
ements for the lowest dimer states that couple to the 0+/^, 
1, and 2 states of SrYb, with the spin-orbit coupling op- 
erator Hso defined within the Breit-Pauli approximation—. 
The spin-orbit coupling matrix elements have been computed 
within the MRCI framework with the MOLPRO code^'*. Di- 
agonalization of the relativistic Hamiltonian gives the spin- 
orbit coupled potential energy curves for the 0+/^, 1 and 2 
states, respectively. Note that all potentials in the Hamilto- 
nian matrices were taken from CCSD(T) and LRCCSD cal- 
culations. Only the diagonal and nondiagonal spin-orbit cou- 
pling matrix elements were obtained with the MRCI method. 
Once the eigenvectors of these matrices are available, one 
can easily get the electric dipole transition moments and the 
nonadiabatic coupling matrix elements between the relativis- 
tic states. In order to mimic the scalar relativistic effects 
some electrons were described by pseudopotentials. For Yb 
we took the ECP28MWB pseudopotential^'', while for Sr 
the ECP28MDF^^ pseudopotential, both from the Stuttgart li- 
brary. For the strontium and ytterbium atoms we used spdfg 
quality basis sets^^'^**, augmented with a set of [2pdfg] dif- 
fuse functions. In addition, this basis set was augmented by 
the set of bond functions consisting of [3s3p2dl f] functions 
placed in the middle of SrYb dimer bond. The full basis of 
the dimer was used in the supermolecule calculations and the 
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Fig. 3 Left: matrix elements of the spin-orbit interaction for the 
a^n, b^L'^, and A'n electronically excited states of SrYb. Right: 
matrix elements of the electric transition dipole moment from the 
'£+ ground electronic state to A'H state (solid curve) and to the 
C'n state (dot-dashed curve), and matrix elements of the 
nonadiabatic angular coupling between the a^Yl and b^Z'^ states of 
SrYb (dashed curve). 

Boys and Bernardi scheme was used to correct for the basis- 
set superposition error . 

Calculations were done for the ground state and 
first fifteen (eight singlet and seven triplet) excited 
states of SrYb. The singlet states correspond to 
the Sr(iD)H-Yb(i5), Sr(ip)H-Yb(i5), STCS)+Yh(4f^5d6s\ 
and Sr('5')H-Yb('P) dissociations, while triplet states to 
Sr(3p)H-Yb(i5), Sr(i5)H-Yb(3p), Sr(3D)H-Yb(i5). The poten- 
tial energies were calculated for twenty interatomic distances 
R ranging from 5 to 50 bohr for each potential curve. The 
ground state potential is presented on Fig. [T] while the poten- 
tial energy curves for the excited states are plotted in Fig. |2] 
The spectroscopic characteristics of all these states are re- 
ported in Table [T] The separated atoms energy for each state 
was set equal to the experimental value. Numeiical values of 
the potentials are available from the authors on request. 

Before going on with the discussion of the potentials let 
us note that the atomic excitation energies obtained from the 
CCSD calculations are accurate. Our predicted position of 
the nonrelativistic state of strontium is 14463 cm^', to be 
compared with the experimental value of 14705 cm^'— de- 
duced from the positions of the states in the P multiplet and the 
Lande rule. For the 'D and 'P states of Sr we obtain less than 
5% difference with the experimental values listed by NIST^". 
Also the life times of the ^Pi and 'Pi states of Sr are accurately 
reproduced. For the ^Pi state we obtained 4.92 ns to be com- 
pared with the expeiimental value of 5.22(3) ns . For the ^Pi 
the theoretical and experimental numbers are 22 and 20 /is^^, 
respectively. Similarly good results were also obtained for the 
ytterbium atom. Such a good agreement between theory and 
experiment for the atoms gives us confidence that the molec- 



Table 1 Spectroscopic characteristics (equilibrium distance, well 
depth, harmonic constant) of the non-relativistic electronic states of 
88gj.l74Yb (jimer. 



State 


Re (bohr) 


De (cm"') 


(Oe(cm ') 


Dissociation 




8.78 


828 


32.8 


Sr('5)+Yb('5) 


A'n 


6.84 


11851 


94.8 


Sr(i£))+Yb(iS) 




7.54 


5201 


63.5 


Sr(iD)+Yb(lS) 


(l)iA 


7.42 


1202 


62.5 


Sr(i£')+Yb(i5) 


(3)iZ+ 


7.91 


2963 


48.5 


Sr(i/')+Yb(i5) 


(2)'n 


7.70 


3112 


61.6 


Sr(i/')+Yb(iS) 


(4)iZ+ 


7.84 


1790 


58.6 


Sr(i/')+Yb(^,|) 


(3)'n 


7.53 


2153 


72.5 


Sr(i/')+Yb(^,|) 


(2)1 A 


7.89 


1175 


40.2 


a^n 


7.02 


6078 


84.7 


Sr(3/')+Yb(i5) 


b3Z+ 


7.84 


4493 


71.3 


Sr(3/')+Yb('5) 


(2)3e+ 
2"'! min. 


7.39 
11.02 


1024 
622 


61.7 
21.0 


Sr(i5)+Yb(3p) 


{2fn 


8.23 


1947 


42.4 


Sr('5)+Yb(V) 


(3)3Z+ 
2"'' min. 


7.45 
9.33 


982 
1077 


92.7 
47.8 


Sr(3£))+Yb(i5) 


(3)3n 


8.04 


1678 


47.9 


Sr(3£))+Yb(iS) 


(1)'A 


7.65 


1422 


50.8 


Sr(^£»)+Yb(iS) 



Table 2 Long-range dispersion coefficients (in a.u.) for ground and 
relevant excited states of the SrYb dimer. 



State 


Ce 


Cg 




2 688 


294 748 


A'n 


3 771 


502 070 




1 265 


509 068 




6 754 


317 656 



ular results will be of similar accuracy, i.e. a few percent off 
from the exact results. 

The ground X'E+ state potential energy curve is presented 
in Fig. [T] It follows from the naive molecular orbital theory 
that the SrYb molecule in the ground state should be consid- 
ered as a van der Waals molecule since the molecular con- 
figuration has an equal number of bonding and antibonding 
electrons. No regular chemical bond is expected, except for 
a weak dispersion attraction and exchange-repulsion. Indeed, 
the ground state potential is weakly bound with the binding 
energy of 828 cm"'. For 7 — it supports A^v = 62 vibra- 
tional levels for the lightest isotope pair, up to A^v = 64 for the 
heaviest isotopes. The changes of the number of bound rovi- 
brational levels and of the position of the last vibrational level 
for different isotopes results in the changes in the scattering 
length from very negative to very positive values. This should 
allow to choose isotopes most suitable for cooling and manip- 
ulation. The equilibrium distance, well depth, harmonic fre- 
quency, and the long-range coefficients of the X'E+ state are 
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reported in Table [T| The permanent dipole moment of SrYb 
in ground electronic state as a function of the interatomic dis- 
tance R is presented in Fig. [T] Except for short interatomic 
distances, the dipole moment is very small. This is not very 
surprising since the two atoms have very similar electroneg- 
ativities and the charge flow from one atom to the other, af- 
ter the formation of the weak van der Waals bond, is very 
small. In fact, similarly as the bonding of the ground state, 
the dipole moment of SrYb should be considered as a disper- 
sion dipole At large interatomic distances it vanishes as 
The vibrationally averaged dipole moment of SrYb 
in ground vibrational state is indeed very small and equal to 
0.023 a.u. 

Potential energy curves of the excited singlet and triplet 
states of SrYb are presented in Fig. |2] An inspection of Fig. 
|2] shows that the potential energy curves for the excited states 
of the SrYb molecule are smooth with well defined minima. 
The potential energy curves of the (2) and (3)^E+ states show 
an avoided crossing, and exhibit a double minimum struc- 
ture. These double minima on the potential energy curves 
are due to the strong nonadiabatic interactions between these 
states. Other potential energy curves do not show any unusual 
features, except for the broad maximum on the potential of 
the (4)'E+ which is most likely due to the interaction with a 
still more excited state not reported in the present work. Ex- 
cept for the shallow double minima of the (2)-^E+ and (3)^1+ 
states, and shallow A states, all other excited states of the SrYb 
molecule are strongly bound with binding energies Dg ranging 
from approximately 1790 cm^' for the (4)'E+ state up to as 
much as 11851 cm^' for the A'H state. 

Let us compare the pattern of the potential energy curves 
of the heteronuclear SrYb molecule with the homonuclear Sr2 
dimer— . In general, molecular orbitals constructed from the 
linear combinations of the Sr(5p)H-Yb(6p) or Sr(4t/)H-Yb(5ii) 
atomic orbitals are expected to have less bonding or antibond- 
ing character than the molecular orbitals constructed from 
the Sr(5p)H-Sr(5p) or ST(4d)+ST(4d) atomic orbitals, because 
large atomic orbital energy differences make combination of 
these orbitals less effective. This would explain why many po- 
tential energy curves of the SrYb dimer are less attractive than 
the corresponding potential energy curves of the Sr2 dimer 
The strongly attractive character of the potential energy curves 
for the first ^E+ and the first ^11 states converging in the long 
range to Sr(^P)H-Yb('5') asymptote could be a result of the sta- 
bilizing effect of the Yh(5d) orbitals for the lowest unoccupied 
orbitals of the a and n symmetry (these molecular are combi- 
nations of the Yh(6p) and Sr(5/?) orbitals), but also of the ap- 
propriate Yh(5d) orbitals, closer in energy to Sr(5 p). Potential 
energy curves for the second ^E+ and second 'll states con- 
verging to Sr(^P)H-Yb('5') are less attractive than the potential 
energy curves for the triplet states, similarly as for the corre- 
sponding states of the Sr2 dimer As for the Sr2 dimer potential 
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Fig. 4 The a-'n, A'n and C'n potential energy curves (solid 

and dashed black curves) in the Hund's case (a) representation that 
are coupled by the spin-orbit interaction and the resulting = 1 
relativistic states (red dotted curves) in the Hund's case (c) 
representation of the SrYb dimer. 



energy curves for the singlet E+ and 11 states converging to the 
Sr(^D)H-Yb(^5) asymptote have a much more attractive char- 
acter than the triplet states converging to the Sr(^D)H-Yb('5) 
asymptote. 

The fl^n, /7^E+, A^n, and C'll excited states essential 
for the photoassociative formation of the ground state SrYb 
molecule proposed in the next section are plotted in Fig. |4] 
The matrix elements of the spin-orbit coupling were calculated 
for the manifolds of coupled a^II, Z7^E+, A'll states, cf. Fig. 
[3] The knowledge of the spin-orbit coupling between a^II, 
/7^E+, A'n, and C'lI states allows us to obtain the relativistic 
(l)O-, (2)0-, (1)0+, (1)1, (2)1,(3)1, (4)1 and (1)2 states by 
diagonalizing the appropriate relativistic Hamiltonian matri- 
ces. The Q.= 1 states are also plotted in Fig. |4] Note that the 
crossing of the Z7^E+ and A'II nonrelativistic states becomes 
an avoided crossing between the (2)1 and (3)1 states. 

Having all the results briefly presented above we are ready 
to discuss the photoassociation process of cold Sr and Yb 
atoms, and look for the prospects of producing ultracold SrYb 
molecules. To conclude this section we would like to empha- 
size that almost all ab initio results were obtained with the 
most advanced size-consistent methods of quantum chemistry: 
CCSD(T) and LRCCSD. In all calculations aU electrons, ex- 
cept for those described by the pseudopotentials, were cor- 
related (42 for ytterbium and 10 for strontium). Only the 
SO coupling matrix elements and the nonadiabatic matrix ele- 
ments were obtained with the MRCI method which is not size 
consistent. Fortunately enough, all of the couplings are impor- 
tant in the region of the curve crossings or at large distances, 
so the effect of the size-inconsistency of MRCI on our results 
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should not be dramatic. 

3 Photoassociation and formation of ground 
state molecules 

Photoassociation is considered for a continuous-wave laser 
that is red-detuned with respect to the intercombination line. 
This transition is dipole-forbidden. However, the a^n state 
correlating to the asymptote of the intercombination line tran- 
sition, cf. Fig.m is coupled by the spin-orbit interaction to two 
singlet states, A^Yl and C'n, that are connected by a dipole- 
allowed transition to the ground electronic state, X ^ E+ . Thus 
an effective transition matrix element is created which can be 
written, to a very good approximation, as 



dso{R) 



{X'-L+\d\C'Yl){C'Yl\Hso\a'n) 



'X^I.+ \d\A^n){A^n\Hso\a^'n) 



(4) 



where H^o is the spin-orbit Hamiltonian in the Breit-Pauli ap- 
proximation^^. The long-range part of dso(R}, dominated by 
the first term in the above expression, is due to the coupling 
with the c'n state, ideally suited for photoassociation. The 
short-range part is due to the coupling with the A^Yl state, 
paving the way toward efficient stabilization of the photoas- 
sociated molecules to the electronic ground state. The a^II 
state, in addition to the spin-orbit coupling with the two sin- 
glet states, is also coupled to the b^T,'^ state correlating to the 
same asymptote, Sr(^P)H-Yb('5). The Hamiltonian describ- 
ing these couplings yielding the Hund's case (c) i2 = 1 states 
reads in the rotating-wave approximation 



H 



(5) 

where H I'^l is the Hamiltonian for nuclear motion in 
the 2S+i|A| electronic state, tf'^'l^l t + r'^'l^l (R) + 

25+1 , 

^trap (R) ^ (1 ^ S„Q)h(OL. The kinetic energy operator is 
given by T = with jj. the reduced mass of SrYb. The 

trapping potential, V^^^^ (R)> is relevant only in the elec- 
tronic ground state for the detunings considered below, even 
for large trapping frequencies. We approximate it by a har- 
monic potential which is well justified for atoms cooled down 
to the lowest trap states. The parameters of the photoassocia- 
tion laser are the frequency, (Ol, and the maximum field ampli- 
tude, E(). The electric transition dipole moments are denoted 
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Fig. 5 Vibrational wave functions of the coupled a-' IT, b^I,^ , A^H, 
and c'n electronic states of SrYb molecule for two binding 
energies corresponding to vibrational quantum numbers v' = — 11 
(left) and v' = — 18 (right) below the dissociation threshold. 

by c/i(R) (X'E+|d|A'n}, t/2(R) = (X'E+|d|C'n), and the 
matrix elements of the spin-orbit coupling are given by 

^i(R) = (fl3n(E = 0,A = ±l)|Hso|/?^E+(E = ±l,A = 0)), 

^2(R) = (fl^'n(E = 0,A = ±l)|Hso|A'n(E-0,A = ±l)), 

^3(R) = (fe^E+(E = ±l,A = 0)|Hso|A'n(E = 0,A = ±l)}, 

^4(R) = (fl^n(E = 0,A = ±l)|Hso|C'n(E = 0,A = ±l)), 

^5(R) = (/7^'E+(E-±l,A-0)|H5o|C'n(E = 0,A = ±l)). 

E and A denote the quantum numbers for the projections of 
electronic spin and orbital angular momenta, S and L, onto 
the internuclear axis. Note that the specific shape of the 
c'n potential energy curve as well as the /^-dependence of 
its spin-orbit coupling and transition dipole matrix elements 
are not important, since the C'n state provides the effec- 
tive dipole coupling only at long range. We have therefore 
approximated the /^-dependence of the couplings with the 
c'n state by their constant asymptotic values in the calcu- 
lations presented below. The Hamiltonian (|5]l has been rep- 
resented on a Fourier grid with adaptive step size (using 

= 1685 grid points and grid mapping parameters /3 = 0.22, 
£^min = 7 • lO^^hartree). 

The key idea of photoassociation using a continuous-wave 
laser is to excite a colliding pair of atoms into a bound level 
of an electronically excited state. ^'^^ For maximum photoas- 
sociation efficiency, the detuning of the laser with respect to 
the atomic asymptote, Sr(-'Pi)H-Yb('5') in our case, is cho- 
sen to coincide with the binding energy of one of the vibra- 
tional levels in the electronically excited state. Fig. |5] shows 
two such levels with binding energies Et, = 5.1cm^' (left) 
and El, = 18.9cm^' (right). Since four electronically ex- 
cited states are coupled by the spin-orbit interaction, the vibra- 
tional wavefunctions have components on all four electroni- 
cally excited states. Note that the norm of the C'H-component 
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Fig. 6 Vibrationally averaged free-to-bound (or 
quasi-bound-to-bound) electric transition dipole moments between 
the lowest trap state of a pair of Sr and Yb atoms colliding in the 
JS^ '£+ ground electronic state in a harmonic trap and two bound 
levels, cf. Figure|5] of electronically excited SrYb dimers as a 
function of the trap frequency. 



of these two vibrational wavefunctions is smaller than 10^^. 
Nevertheless, this is enough, similar to the photoassociation 
of the strontium dimers near an intercombination line,— to 
provide the transition dipole for the free-to-bound (or quasi- 
bound-to-bound, due to the trapping potential) excitation. The 
vibrational level with binding energy Ei, = 5.1 cm^' is pre- 



dominantly of triplet character (with 56% of its norm resid- 
ing on the fl^n state, 32% on the b^Y.^ state and just 11% 
on the A'n state), while the vibrational level with binding 
energy Ef, — 18.9cm^' shows a truly mixed character (55% 
triplet vs 45% singlet). The fact that multiple classical turn- 
ing points are clearly visible in the vibrational wavefunction 
with Ef, = 18.9cm^^ reflects the resonant nature of the spin- 
orbit coupling of this level: the coinciding energy of the levels 
in the coupled vibrational ladders leads to a resonant beating 
between the different components of the coupled wavefunc- 
tions.^*^ Such a structure of the vibrational wavefunctions was 
shown to be ideally suited for efficient stabilization of the pho- 
toassociated molecules into deeply bound levels in the ground 
electronic state. — ,"Zi, 

The Condon radius for photoassociation coincides with the 
classical outer turning point, i.e., roughly speaking with the 
outermost peak of the vibrational wavefunctions as shown in 
Fig. |5] Since the pair density of the atoms colliding in their 
electronic ground state decreases with decreasing interatomic 
distance, photoassociation is more efficient for small detun- 
ing. This is reflected by the larger values of the black com- 
pared to the red curve in Fig. |6]which shows the free-to-bound 
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Fig. 7 Vibrationally averaged bound-to-bound electric transition 
dipole moments between the vibrational levels of the coupled 
electronically excited states that are shown in Fig.|5]and all 
vibrational levels of theX'Z+ ground electronic state. 



(quasi-bound-to-bound) transition matrix elements for the two 
vibrational wavefunctions depicted in Fig. |5]as a function of 
the trapping frequency of the optical lattice. The second ob- 
servation to be drawn from Fig. |6]is the almost linear scaling 
of the transition matrix elements, and hence the photoassoci- 
ation probability, with the trap frequency. That is, enhancing 
the trap frequency from 50kHz, which has been employed for 
photoassociation of Sr2^^, to 500 kHz, which is within cur- 
rent experimental feasibility, will increase the number of pho- 
toassociated molecules by a half an order of magnitude. This 
confinement effect is easily understood in terms of the larger 
compression of the quasi-bound atom pairs in a tighter optical 
trap. 

In view of the formation of deeply bound molecules in their 
electronic ground state, it might be advantageous to choose 
the larger detuning of 18.9cm^' despite the photoassociation 
probability being smaller by about a factor of 5.9 compared 
to a detuning of 5.1 cm^' for all trap frequencies. This be- 
comes evident by inspecting Fig. [T] which displays the bound- 
to-bound transition matrix elements between the two electron- 
ically excited vibrational wavefunctions with Ef, = 5.1 cm^' 
and Ef, = 18.9cm^^ and all bound levels of the X'E+ elec- 
tronic ground state. These transition matrix elements govern 
the branching ratios for spontaneous decay of the photoasso- 
ciated molecules. Note that for v' = — 11 and v' ~ —18, the 
electronically excited molecules will decay into bound levels 
of the electronic ground state with a probability of about 24%. 
This decay to a large extent into bound levels is a hallmark 
of photoassociation near an intercombination line.^^ It is in 
contrast to photoassociation using a dipole-allowed transition 
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where the probabihty for dissociative decay is often several or- 
ders of magnitude larger than that for stabilization into bound 
ground state levels. ^'^ While the excited state vibrational level 
with El, = 5.1 cm^' has its largest transition dipole matrix el- 
ements with the last bound levels of the X'E+ ground elec- 
tronic state that are only weakly bound, a striking difference 
is observed for the excited state vibrational wavefunction with 
iii = 1 8 .9 cm^ ' . The strong singlet-triplet mixing of this level, 
in particular the pronounced peak near the outer classical turn- 
ing point of the A^Yl state, cf. Fig. |5] leads to significantly 
stronger transition dipole matrix elements with deeply bound 
levels of the X^E+ ground electronic state for v' = — 18 com- 
pared to v' = — 11, the one with v" = 1 being the largest. Of 
course, the transition dipole matrix elements govern not only 
the spontaneous decay of the photoassociated molecules but 
also stabilization via stimulated emission. Due to the compar- 
atively long lifetime of photoassociated molecules, estimated 
to be of the order of 15 /is, stabilization into a selected single 
vibrational level of the electronic ground state can be achieved 
by stimulated emission using a second continuous-wave laser 

Before outlining how a prospective experiment forming 
SrYb molecules in their vibronic ground state based on our 
results could proceed, it is natural to ask whether the accuracy 
of the calculations is sufficient for such a prediction. In par- 
ticular, how sensitively do our results for the binding energies 
and structure of the vibrational levels as well as for the transi- 
tion matrix elements depend on the accuracy of the electronic 
structure calculations? The binding energies depend mostly 
on the quality of the potential energy curves, where the er- 
ror is estimated to be a few percent, and to some extent, for 
the spin-orbit coupled excited states, on the accuracy of the 
spin-orbit interaction matrix elements (error of a few percent). 
Nevertheless, the uncertainty of our potential energy curves 
is smaller than the range of reduced masses, as illustrated in 
Fig. [U Therefore photoassociation with subsequent stabiliza- 
tion to a low-lying vibrational level should work for all isotope 
pairs since levels with strong perturbations due to the spin- 
orbit interaction are always present in the relevant range of 
binding energies, respectively, detunings, cf. Fig.|8] 

In fact, the exact position and the character of the excited 
state vibrational level, strongly perturbed such as the one 
with Et, = 18.9cm^' or more regular such as that with E/, = 
5.1 cm^' in Fig. |5] can be determined experimentally. ^^'^^ A 
possible spectroscopic signature of the character of the vibra- 
tional wavefunctions is the dependence of the rotational con- 
stants, {v'\ ^^'^T \ v'), on the binding energy of the correspond- 
ing levels. This is shown in Fig.|8]for different isotope com- 
binations of strontium and ytterbium. The rotational constants 
of those levels that are predominantly of triplet character lie on 
a smooth curve, while those that are mixed deviate from this 
curve. This behavior is easily understood as follows: Without 
the coupling due to the spin-orbit interaction, the rotational 
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Fig. 8 Rotational constants of the vibrational levels of the coupled 
a^n, A^n, and C'n electronically excited states of the SrYb 

molecule for different isotope pairs. The isotope **^Sr'^^Yb was 
employed in the calculations shown in Figs. l5l7l andl9l 



constants of the a^Yl, b^T.^ and A'll states would each lie 
on a smooth curve. The strongly mixed levels 'belong' to all 
three curves at the same time. Correspondingly, the value of 
their rotational constant lies somewhere inbetween the smooth 
curves of the regular levels. The lower peaks at small binding 
energies in Fig.[8]indicate mixing mostly between the a^II and 
/7^E+ states, while the high peaks at larger binding energies 
reflect a strong singlet-triplet mixing. Spectroscopic determi- 
nation of the rotational constants thus allows for identifying 
those excited state levels that show the strongest singlet-triplet 
mixing and are best suited to the formation of ground state 
molecules. Spectroscopy is also needed to refine the value for 
the transition frequency of the stabilization laser The binding 
energies of vibrational levels of the X'E+ ground electronic 
state come with an error of 5 %, i.e., ± 50 cm^\ defining the 
window for spectroscopic search. 

Note that our model, Eq. (|5]l, does not account for angular 
couplings, i.e., couplings of the £2=1 states with £2 = 0^ 
and Q. — 2. When including non-adiabatic angular couplings, 
we found the components of the vibrational wavefunctions on 
the newly coupled surfaces to account for less than 0.001% 
of the population. The changes in the binding energy of the 
vibrational levels turned out to be less than 10^''cm^', well 
within the error of the electronic structure calculations. This 
negligible effect of the angular (Coriolis-type) couplings for 
SrYb is not surprising due to its large reduced mass which 
enters inversely all coupling matrix elements. 

Combining all results shown above and assuming that the 
relevant spectroscopic data has been confirmed or adjusted ex- 
perimentally, we suggest the following scheme for photoasso- 
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Fig. 9 Proposed scheme for the formation of ground state SrYb 
molecules via photoassociation near the intercombination line 
transition with detuning Ag,^ = 18.9cm"' (Vtrap = 100 MHz). 

ciation of SrYb dimers followed by stabilization via stimu- 
lated emission (see Fig.|9]i: 

1 . A large trapping frequency of the optical lattice is chosen 
to optimally compress the pair density of strontium and 
ytterbium atoms prior to photoassociation. 

2. A photoassociation laser with frequency Oi ~ 690 nm, 
red-detuned from the intercombination line transition and 
resonant with an electronically excited vibrational level, 
v', of strongly mixed singlet-tripled character, is applied 
for a few jtis. The duration of the photoassociation laser 
(about 5 /is roughly is an upper bound) is a compromise 
between saturating photoassociation and avoiding spon- 
taneous emission losses (lifetime of about 15 /is) while 
the laser is on. 

3. As the photoassociation laser is switched off, the stabi- 
lization laser is switched on. Due to the strong bound-to- 
bound transition matrix elements, saturation of the tran- 
sition is expected already for shorter pulses (< 1 /is). The 
frequency of the stabilization laser, ©2 ~ 655 nm, is cho- 
sen to be resonant with the transition from the electroni- 
cally excited level, v', to the first excited vibrational level 
of the X'E+ electronic ground state, v" = 1. 

4. Before repeating steps 2 and 3, both photoassociation 
and stabilization lasers remain turned off for a hold pe- 
riod in which the X'E+(v" = 1) molecules decay to 
the vibronic ground state, X'E+(v" = 0). This en- 
sures that the molecules created in the electronic ground 
state by the first sequence of the photoassociation and 



stabilization steps are not re-excited in a following se- 
quence. The formed molecules can then be accumulated 

mX^T.+ {v" = 0). 

Note that this scheme does not require phase coherence be- 
tween the two pulses. Step 4 needs to involve a dissipative el- 
ement in order to ensure the unidirectionality of the molecule 
formation scheme. Dissipation can be provided by infrared 
spontaneous emission due to the permanent dipole moment of 
the heteronuclear dimers. However, this timescale is estimated 
to be of the order of 5 s, much too slow to be efficient for ac- 
cumulation of ground state molecules. A second possibility is 
due to collisional decay. For the decay to occur within 1 ms, 
a density of lO'^cm"^ is required. Note that the density was 
3 • 10'^ cm"^* in the experiment photoassociating Sr2 in an op- 
tical lattice with trapping frequency 50kHz. Increasing the 
trap frequency will further increase the density such that hold 
times in the sub-ms regime are within experimental reach. 

One might wonder whether the comparatively long hold 
times can be avoided by using Stimulated Raman Adiabatic 
Passage (STIRAP)^^ for the photoassociation (pump) and sta- 
bilization (Stokes) pulses. ^^'^^ In order to overcome the prob- 
lem of unidirectionality that occurs in repeating the photoasso- 
ciation and stabilization steps many times, the whole ensemble 
of atom pairs in the trap needs to be addressed within a single 
STIRAP sweep or within a single sequence of phase-locked 
STIRAP pulse pairs. Note that the Stokes/stabilization pulse 
should be tuned to the v' — > v" = transition in this case. The 
feasibility of STIRAP-formation of ground state molecules 
depends on isolating the initial state sufficiently from the scat- 
tering continuum. A possibility to achieve this that was dis- 
cussed theoretically consists in utilizing the presence of a Fes- 
hbach resonance. ^^'^^ If no resonance is present, i.e., in an 
unstructured scattering continuum, STIRAP fails. In a se- 
ries of ground-breaking experiments, STIRAP transfer to the 
ground state was therefore preceded by Feshbach-associating 
the molecules.ii^"— An alternative way to isolate the initial 
state for STIRAP from the scattering continuum that does not 
rely on Feshbach resonances (which are absent for the even 
isotope species of Sr and Yb) is given by strong confine- 
ment in a deep optical lattice. An estimate of the required 
trap frequency is given in terms of the binding energy of the 
Feshbach molecules that were STIRAP-transferred to the vi- 
bronic ground state. It was for example about 230kHz for 
KRb molecules. Hence a deep optical lattice with trapping 
frequency of the order of a few hundred kHz should be suffi- 
cient to enable STIRAP-formation of ground state molecules. 
However, in order to be adiabatic with respect to the vibra- 
tional motion in the trap with periods of the order of about 
1 /IS, the duration of the photoassociation pulse needs to rather 
long, at least of the order of 10 /is. The challenge might be to 
maintain phase coherence between the photoassociation pulse 
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and the stabilization pulse over such timescales. For a train of 
phase-locked STIRAP-pulse pairs, the requirement of dura- 
tions of the order of 10 /xs or larger applies to the length of the 
sequence of pulse pairs. The minimum Rabi frequencies to 
enforce adiabatic following are = 159 kHz for a 10 /xs-pulse 
or £2 = 15.9kHz for a 100/is-pulse. As a further prerequi- 
site, all or at least most atom pairs should reside in the lowest 
trap state, Vtrap = 0. Then steps 2-4 above might be replaced, 
provided the trapping frequency is sufficiently large, by 

2.' a single STIRAP-sweep^- forming ground state 
molecules with /xs-pulses where the stabilization laser, 
tuned on resonance with the v' ^> v" = transition 
((02 ~ 654nm), precedes the photoassociation laser, 
tuned on resonance with the Vtrap = ^ v' transition 
(tO2«690nm); 

2." or, a train of short, phase-locked STIRAP pulse pairs with 
correctly adjusted pulse amplitudes.— 

To convert the Rabi frequencies to field amplitudes, note 
that the transition matrix elements are 5 10-^ for the pump 
pulse (assuming a trap frequency of 300kHz) and 3 • 10^^ for 
the Stokes pulse. Phase coherence needs to be maintained 
throughout the single STIRAP-sweep or sequence of STIRAP 
pulse pairs. 

4 Summary 

Based on a first principles study, we predict the photoassocia- 
tive formation of SrYb molecules in their electronic ground 
state using transitions near an intercombination line. The po- 
tential energy curves, non-adiabatic angular coupling and spin 
orbit interaction matrix elements as well as electric dipole 
transition matrix elements of the SrYb molecule were calcu- 
lated with state-of-the-art ab initio methods, using the coupled 
cluster and multireference configuration interaction frame- 
works. Assuming that the accuracy of the calculations for the 
SrYb molecule is about the same as for the isolated Sr and Yb 
atoms at the same level of the theory, we estimate the accuracy 
of the electronic structure data to 5%. However, the crucial 
point for the proposed photoassociation scheme is the exis- 
tence and position of the intersection of the potential energy 
curves corresponding to b^T, and A^H states. By contrast to the 
binding energies of the vibrational levels, the position of this 
intersection does not depend very much on the overall quality 
of the computed potential energy curves curves. The correct 
structure of the crossings between the potential curves of the 
a^Yl, b^l, and A'n states is reproduced using even relatively 
crude computational methods of quantum chemistry which do 
not account for dynamic correlations such as the multiconfig- 
uration self-consistent field (MCSCF) method employed here. 



The spin-orbit coupled a^Yl, fe^E+, A^H, and C'n electron- 
ically excited states are essential for the photoassociation. A 
pair of colliding Sr and Yb atoms is excited into the triplet 
states with the dipole coupling for the photoassociation (sta- 
bilization) provided by the C^Yl (A^H) state. The formation 
of SrYb molecules in their electronic ground state proceeds 
via photoassociation into a weakly bound level of the coupled 
electronically excited states (Oi w 690 nm) followed either by 
spontaneous or stimulated emission. 

If photoassociation is followed by spontaneous emission, 
about 24% of the photoassociated molecules will decay into 
bound levels of the ground electronic state, roughly indepen- 
dent of the detuning of the photoassociation laser However, 
which ground state rovibrational levels are populated by spon- 
taneous emission depends strongly on the detuning of the pho- 
toassociation laser While most detunings will lead to decay 
into the last bound levels of the ground electronic states, cer- 
tain detunings populate excited state levels with strong spin- 
orbit mixing. The strongly resonant structure of the wave- 
functions allows for decay into low-lying vibrational levels. 
This might be the starting point for vibrational cooling if 
molecules in their vibronic ground state are desired. 

Alternatively, the long lifetime of the photoassociated 
molecules, of the order of 15jUs, allows for stabilization to 
the electronic ground state via stimulated emission, by a se- 
quence of photoassociation and stabilization laser pulses of 
jUs duration. Two schemes are conceivable: (i) A repeated 
cycle of photoassociation and stabilization pulses is applied 
with X^E+(v" = 1) as the target level. The duration of the 
pulses should be of the order of 1 /is. In order to accumu- 
late molecules in X^E+(v = 0), a hold period whose duration 
depends on the density of atoms is required for collisional de- 
cay from V = 1 to V = 0. For deep optical lattices, hold pe- 
riods in the sub-ms regime can be reached, (ii) The vibronic 
ground state, X'E+(v — 0), is targeted directly by a counter- 
intuitive sequence of photoassociation and stabilization pulses 
(STIRAP), either using two long pulses or a train of phase- 
locked pulse pairs. The timescale for the pulses is deter- 
mined by the requirement to be adiabatic with respect to the 
motion in the optical lattice. The largest trapping frequencies 
feasible to date imply pulse durations at least of the order of 
10/is. Phase coherence between the pulses needs to be main- 
tained over this timescale. Note that STIRAP fails if applied 
to an unstructured scattering continuum of colliding atoms. 
A possibility to circumvent this is given by preselecting the 
initial state for STIRAP with the help of a (Feshbach) reso- 
nance. ^^"^^ Our variant of the scheme is different since STI- 
RAP is enabled by the presence of a deep trap. 

Before either of the above discussed molecule formation 
schemes can be implemented experimentally, our theoretical 
data needs to be corroborated by spectroscopy. In particular, 
our binding energies come with an error of a few percent, im- 
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plying a corresponding uncertainty in the transition frequen- 
cies. Moreover, the exact position of strongly spin-orbit mixed 
excited state wavefunctions needs to be confirmed by mea- 
suring the excited state level spacings or rotational constants. 
However, despite the relatively large uncertainties in the ener- 
gies of the rovibrational levels important for the proposed pho- 
toassociation scheme, our ab initio methods correctly locate 
the crossing of the singlet and triplet potential energy curves. 
This is the key ingredient for the efficient production of ground 
state SrYb molecules that we are predicting with our study. 
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